Railway switch transport model 
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We propose a simple model of coupled heat and particle transport based on a zero-dimensional 
classical deterministic dynamics which is reminiscent of a railway switch whose action is only a 
function of the particle's energy. It is shown that already in the minimal three-terminal model, 
where the second terminal is considered as a probe with zero net particle and heat currents, one can 
find extremely asymmetric Onsager matrices as a consequence of time-reversal symmetry breaking 
of the model. This minimalistic transport model provides a better understanding of thermoelectric 
heat engines in the presence of time-reversal symmetry breaking. 
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Introduction. - Minimalistic mathematical models often 
provide key paradigms in theoretical physics. In the the- 
ory of coherent quantum transport [1 of non-interacting 
electron systems, the conductances can be elegantly for- 
mulated solely in terms of transmission and reflection 
amplitudes . Very similar expressions for conductances 
can be written also in the realm of classical physics, if the 
particle's dynamics inside the systems is deterministic 
and conservative (Hamiltonian) (see e.g. [3]). In multi- 
terminal transport theory, probe terminals |3] are conve- 
niently used as minimalistic models of inelastic scatter- 
ing. A probe is a terminal whose temperature and chem- 
ical potential is chosen self-consistently so that there is 
no average flux of particles and energy between the probe 
itself and the system under examination. The advantage 
of such approach lies in its simplicity and independence 
of microscopic details of inelastic processes. Probe termi- 
nals have been widely used in the literature and proved 
to be useful to unveil nontrivial aspects of phase-breaking 
processes [I] , heat transport and rectification [5HTT] , and 
thermoelectric transport [T2HT8] . 

In this Letter we propose and study a simple multi- 
terminal classical transport model, where the wires con- 
nected to the system are one-dimensional, i.e. they have 
just one momentum state per each value of the energy, 
and where the deterministic dynamics of the system, at 
fixed energy, is simply a permutation (rewiring) among 
the terminals. In particular, we will focus on three- 
terminal models, depicted in Fig. [l] where the second 
terminal is considered as a probe. We discuss the cou- 
pled heat and particle deterministic transport between 
the remaining wires within the linear response limit. In 
this limit the transport between the first and the third 
wire is described by the reduced Onsager matrix. As our 
energy-dependent rewiring is in general not self-inverse 
we are particularly interested in the precise conditions 



for the reduced Onsager matrix to be asymmetric and 
discuss how to maximize the asymmetry. 
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FIG. 1: Schematic drawing of the three terminal transport 
model, with jp^i and jq,i denoting the particle and heat cur- 
rents from terminal i. 

Classical coupled transport formalism.- Let us consider 
a generic A'^-terminal, non-interacting classical transport 
model. Each wire labelled by i is connected to a thermo- 
chemical bath at reciprocal temperature Pi = fi^5(3i and 
chemical potential Hi = jj, -\- SjjLi {i = 1, ...,N), where /3 
and fj, are reference reciprocal temperature and chemical 
potential, respectively. The S/3i and Sjii are considered as 
gradient fields w.r.t. reference values. Such gradients are 
in the linear response regime small in magnitude. The 
particles are effused from the wires into the junction, with 
the injections rates 

= (1) 

Pi 

where 7 denotes the injection rate at the reference values 
/3 and /i. Without loss of generality, we set j — (3 — 1, 
and = 0. At inverse temperature /3, the particles 
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energies are distributed according to Boltzmann's for- 
mula p{E) = [ie~^^ . In the stationary, linear response 
regime, we have in the iih wire a particle current jp^i 
and a heat current jq^i proportional to gradients 5[3j 
and —(35^j = —Sjij [T5]. To state this more precisely, 
we introduce the 2- vector of currents Ji = {jq^i,jp,i) 
in the ith terminal, and the 2-vector of gradient fields 
Xj = (SPjj—SjjLj) in the jth terminal, and 2x2 blocks 
of the Onsager matrix 



connecting the two of them: 

N 

The matrix elements of Lij are defined by 



dEe 



-E~ 



■\E){l,E,E^ 



(2) 



(3) 



(4) 



with fj j = 5i_j — Ti^j, and the on-shell transmis- 
sion functions Ti^j satisfying the probability conserva- 
tion = 1 ^^'^ the sum rule J2i'''i^3i-^) = 
Ti^^j (E) , ensuring that the currents vanish when all 
the potentials and temperatures are equal. The conserva- 
tion of the net currents translates to a 2-vector condition 
X^ili '^i — 0- We note that due to conservative and non- 
interacting nature of our model, the Onsager matrix is 
always, as explicitly written in Eq. ([2]), block symmetric. 
This property is a consequence of conservation of total 
probability and would translate also in any quantum ex- 
tension of our model, being in that case a consequence of 
the unitarity of the S-matrix . 

The railway switch model.- We consider a determinis- 
tic transport model in which the outgoing wire for a par- 
ticle is uniquely determined by the incoming wire and 
the particle's energy. Such model is not bound to be 
symmetric with respect to time reversal, and therefore 
allows us to systematically study the effects of time re- 
versal symmetry breaking. The model is fully specified by 
transmission functions Ti<_ j (E) for which only the values 
zero and one are allowed: Ti^j{E) = 1 if particles in- 
jected from terminal j with energy E go to terminal i, 
Ti^j{E) = if such particles go to a terminal different 
from i. The above deterministic on-shell transmission 
functions can be described using permutation matrices 
Pk for k E I = {1, . . . , iV!} corresponding to the permu- 
tation group Sn as 



(E) = [P^^E)h 



i,j = 1,...,N. 



(5) 



where ip{E) : M ^ / is a piecewise constant function con- 
trolling the (energy-dependent) switching between the 
permutations. The function ipiE) for n permutation 



switches is completely specified in terms of a sequence 
of n -I- 1 integers tt = {po,Pi, ■ ■ ■ ,Pn), Pi G I, and a se- 



quence of n threshold energies e = {Ei 
the switches occur: 



, En), at which 



^P{E)=p,El for Ee[E,,E,+{\ 



(6) 



with i = 0, . . . , n, i?o 
at energy E^ (« = 1, 



and En+i = oo. That is to say, 
we switch from permutation 



Pp-_^ to Pp-. A realization of the model is defined by a 
pair (tt, e). 

From now on we focus on the (N — 3)-terminal case 
(Fig.[l]), so that the permutation matrices Pk G {0, Ij^^s 
read 



Pi 



Ph 



P^ 



1 
1 
1 

1 0' 

1 
1 

1 

1 
1 



V 



P2 



Pa 



1 
1 
1 

1 0' 

1 

1 

1 

1 

1 



X7- 



To the right of each permutation we show the rewiring 
described by the matrix. Following the figures we see 
that each permutation has a very simple interpretation, 
e.g., Pi represents the refiection of the incoming particles 
to the original terminals, in P2 particles from terminal 1 
are reflected, particles from terminal 2 go to terminal 3 
and particles from terminal 3 go to terminal 2, etc. 

Our model has interesting symmetries. By introducing 
the super-operator of time inversion I'X — X"'" (with X^ 
the transpose of X), we can see that 

fPk^Pk for A: €{1,2,3,6} and fp^^P^. (7) 

This implies that P4 and P5 are the only time-asymmetric 
permutations, with P4 = P^^ and P5 = P^^ , while for 
all other cases we have P^ — P^^ ■ By considering the 
second terminal to act as a probe, the model is invari- 
ant on swapping the first and the third terminal. These 
operation can be performed by the swap super-operator 
SX = PqXPq acting as 



SPi = Pi 



SPo = P. , SPi = p. 



SPe = Pe . (8) 



Thermoelectric transport with a probe terminal.- To il- 
lustrate the railway switch model in a concrete example, 
we discuss thermoelectric (or thermochemical) transport, 
when the second terminal acts as a probe, i.e., J2 = 0. 
We choose to measure gradients w.r.t. the third wire 
by setting X^, = 0. We end up relating the remaining 
gradient fields 



X. 



L22 L21X1, 



(9) 
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and connecting the currents between the first and the 
third wire, 

■h = L^cdXi, (10) 
in terms of a reduced Onsager matrix Lrcd, defined as 

^rcd = Lii — L12L22 L21 ~ 



hi h2 
hi h2 



(11) 



The reduced Onsager matrix admits a nice analytic de- 
scription of matrix element h.j- Let us first introduce the 
determinant D = dot L22 of the Onsager matrix L22 for 
the transmission from the probe, written as 



D 



1 



Yle-''^f22iE,)dE, 



{Ei-E2f (12) 



and define the integration measure 



d/i 



n 



e-^^AE, 



T2i{E^)ti2{E2)t22{Ez) (13) 



over the domain {Ei, E2, E3) e V = R^. (To simphfy 
notation we omit the arguments (Ei, £^2,-^3) in the mea- 
sure dfi.) Then we can write the matrix elements of Li-od^ 



Kii - ^ J^dfiEiE2{Ei - E3){E2 - E3)IU) 



hi 



h2 = T11-- dfi {El ~ E^)[E2 - ^3), (15) 
h2 = Qu^ ^ J^df,fiEi,E2,E3), 



21 



Qii-^J^dfif{E2,Ei,E3), 



(16) 
(17) 



with f{Ei,E2,E3) = E2EI EiEl - {El + E2)E2E3. It 
is convenient to introduce the quantities 

2/ := f{Ei,E2,E3) + f{E2,Ei,E3) 

= {El + E2){Ei ~ E3){E2 - E3) , (18) 
2A/ := f{Ei,E2,E3)- f{E2,Ei,E3) 

= ^e.,,,fei?,i?L (19) 

with j.fc being the Levi-Civita totally asymmetric ten- 
sor. Note that L^cd can be decomposed into the sum of 
Lii and a matrix defining the communication between 
the wires. 

In this work we are mainly interested in the asymmetry 
measure x of the reduced Onsager matrix, defined as 



hi 

h2 



(20) 



This parameter has been discussed for quantum dots with 
broken time-reversal symmetry |15[ I21j . If the trans- 
port is time-reversal symmetric, then J d/zA/ = and 



consequently x = 1. Large asymmetries are desirable 
since they could in principle lead to high thermoelec- 
tric efficiencies. Indeed, by considering the reduced sys- 
tem as a heat engine with steady state power generation 
P — —jp.iS^i, dissipated heat current Q = jq^i, and 
efficiency rj = P/Q, the maximum efficiency reads [22] 



Vn 



V 



1 



(21) 



where rye — \5Pi\ (for /? = 1) is the Carnot efficiency 
and y = je^jf is the figure of merit. Note that x acts 
as a multiplier to the efficiency and so its maximization 
(at fixed y) is desirable for increasing the efficiency. The 
thermodynamic bounds on x, y, and 77 are discussed in 
Ref. [22 . 

In the following wc discuss the asymmetry measure x 
and the figure of merit y of the reduced Onsager matrix in 
different realizations of the railway switch model model. 
At fixed threshold energies e, by considering the sequence 
of inverted permutations we obtain the same figure of 
merit and reciprocal measure of asymmetry |23j : 



y{fTT) = y(7r) 



x{Ttt) = x{'k) 



(22) 



while exchanging the first and the third wire does not 
affect the two quantities: 



y{Sii) = y{-K) 



x{Stt) = x{Tr) 



(23) 



If a realization is symmetric w.r.t. to tt — Ttt, then is has 
a symmetric Onsager matrix and so a; = 1. We observed 
numerically that the Onsager matrix of our model has 
all elements positive, h,j > 0, so that the asymmetry 
measure x is always non-negative, x > 0. 

The number of all n-switch cases is Nn = 6-5". Let 
Sn denote the number of cases for which we have a; = 1. 
As shown in Fig. [2] 



as n 



00 , 



(24) 



with a « 3 and C ~ 24. We see that Sn is asymptotically 
approximately six times larger than the number of cases 
composed of only symmetric permutations, equal to Sn — 
4-3". The number of cases An = Nn — Sn with an 
asymmetric Onsager matrix increases exponentially, as 
0(5"). Hence, most of the cases lead to a; 7^ 1. 

In the following, we discuss the properties of our model 
for different number n of switches, thus increasing with 
n the complexity of the model in a controllable manner. 
All one-switch (n = 1) cases have x = I. The asym- 
metry X ^ 1 is possible only in the cases with two or 
more switches. In the two-switch (n — 2) cases only non- 
repeated combinations of permutations Pk with indices in 
k g {2,3,4,5} produce asymmetry, but the asymmetry 
parameter x is always limited to a finite interval, namely 
it is always strictly larger than zero and finite. In the 
three-switch {n = 3) cases the asymmetry parameter can. 
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FIG. 2: The number of cases Sn with symmetric Onsager 
matrix as a function of the number n of switches and best 
fitted asymptotic dependence 24-3". 

for specific sequences of permutations, become arbitrar- 
ily large. This fact is possible for the following sequences 
of permutations: 

^ = {(2,3,1,4), (3,2,1,5), (4,2,1,3), (25) 
(4,2,1,5), (5,3,1,2), (5,3,1,4)}. 

The obtained set is invariant w.r.t. S, swapping the first 
and the third channels. Therefore we can limit ourselves 
to consider a desymmetrized set 

y = {(2,3,1,4), (4,2,1,5), (4,2,1,3)}. (26) 

Equally interesting are the cases in which x limits to 
0, which are obtained by time-inverting (tt — > Ttt) the 
cases in V. In all these cases we can tend to the maximal 
asymmetry, provided the switch energy thresholds are 
chosen with certain asymptotic scalings. By expressing 
the threshold energies S^+i = Ei + AE^ in terms of the 
gaps AEi > 0, for the cases of V, a local maximum Xmax 
of X at fixed AE2 follows a curve scaling in the limit 
AE2 — >■ 00 as 

Eo-e-"^^\ A£;i X e-^^■^^ (27) 

and asymmetry x diverges along this curve according to 

x„.ax >^ e("-«^^= . (28) 

For the first two cases in V we find (a, (3) — 
(0.488,0.316), for the last one (a,/3) = (0.842,0.536). 
We note that three-switch cases of V are obtained from 
two-switch asymmetric cases by inserting permutation Pi 
into the third position. We maximize the asymmetry by 
increasing the energy interval controlled by Pi permuta- 
tion and decrease all other intervals. 

In Fig. [3] we show how the figure of merit y is related 
to the asymmetry measure x in two examples, one from 
TT e y and the other corresponding to its time-inversion 
TT — > Ttt. We see that increasing the asymmetry x results 
on average in decreasing the figure of merit y, in such 




2 4 6 8 10 



FIG. 3: The relation between x and y for random choices of 
threshold energies in the case (2, 3, 1, 4) (red dots) and in the 
time inverted case of (4, 2, 1, 3), that is, in the case (5, 2, 1, 3) 
(green dots). The black curve gives the analytic boundary 
y = 4x/(x — 1)^. 

a way that we do not improve the efficiency ry. Note 
that, for a; > 1, we are far from the Carnot efficiency, 
which is achieved on the curve h{x) — 4a;/(x — 1)^ [25] 
(black curve in Fig. [s]). Further increasing the number 
of switches the diversity of cases increases beyond the 
ability of a detailed case-by-case study. It remains an 
open problem whether one could with our model and for a 
large number of switches, approach the Carnot efficiency 
at asymmetries x 3> 1. 

Generalizations - By replacing permutation matri- 
ces Pi G 53 with matrices corresponding to the per- 
mutation group of arbitrary degree N , our model de- 
scribes deterministic transport in an iV-terminal junction 
of one-dimensional wires. Seeking optimized transport 
in our model therefore corresponds to discrete optimiza- 
tion problems on permutation groups. As a principal, 
but completely impractical generalization we note that 
we can phrase any deterministic scattering dynamics in 
terms of our switch model. Finally, our model could be 
easily extended to investigate nonlinear regimes where 
breaking of time-reversibility has nontrivial effects on the 
transport [24], and reformulated in a quantum mechani- 
cal context. 

Conclusions - We have presented a minimalistic clas- 
sical finite-state deterministic transport model which 
allows to systematically study the questions of cou- 
pled particle and heat transport and thermoelec- 
tric/thermochemical efficiency. We have analyzed in 
particular the three-terminal model, where one terminal 
serves as the probe, and analyzed in detail the conditions 
under which one can maximize the asymmetry of the re- 
duced Onsager matrix in relation to time-reversal break- 
ing in the model. We expect that our model may serve 
as a simple benchmark for mesoscopic coupled transport 
studies. M.H. and T.P acknowledge supported by the 
grant Pl-0044 of Slovenian Research Agency, and G.B. 
and G.C. by the MIUR-PRIN 2008 and by Regione Lom- 
bardia. 
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